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Abstract 

We elucidate constraints imposed by confinement and dynamical chiral symmetry 
breaking on the infrared behaviour of the dressed-quark and -gluon propagators, and 
dressed-quark-gluon vertex. In covariant gauges the dressing of the gluon propagator 
is completely specified by V(k 2 ) := 1/[1 + IT(A; 2 )], where H(k 2 ) is the vacuum polari- 
sation. In the absence of particle-like singularities in the dressed-quark-gluon vertex, 
extant proposals for the dressed-gluon propagator that manifest V(k 2 = 0) = and 
max(7-'(A; 2 )) ~ 10 neither confine quarks nor break chiral symmetry dynamically. 
This class includes all existing estimates of V{k 2 ) via numerical simulations. 
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Strong interaction phenomena are characterised by dynamical chiral symme- 
try breaking (DCSB) and colour confinement. At low energy, DCSB is the 
more important; for example, in its absence the tt- and p-mesons would be 
nearly degenerate and at the simplest observational level that would lead to a 
markedly different line of nuclear stability. These phenomena can be related to 
the infrared behaviour of elementary Schwinger functions in QCD and herein 
we elucidate some constraints they place on this behaviour. 

As described pedagogically in Ref. [1], DCSB can be studied using the QCD 
"gap equation"; i.e., the Dyson-Schwinger equation (DSE) for the renormalised 
dressed-quark propagator (connected, 2-point, dressed-quark Schwinger func- 
tion), S(p): 
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= Z 2 (i7 -p + m bm ) + y g 2 D fll/ {p- q)l^S{q)T v {q,p) , (2) 



where D^k) is the renormalised dressed-gluon propagator, r M (g;p) is the 
renormalised dressed-quark-gluon vertex, and J* g A := j A d 4 q/(27r) 4 represents 
mnemonically a translationally-invariant regularisation of the integral, with 
A the regularisation mass-scale. The final stage in any calculation is to take 
the limit A — > oo. In (1), Z\ and Z 2 are the renormalisation constants for the 
quark-gluon vertex and quark wave function, and m^, m is the current-quark 
bare mass: the chiral limit is obtained with mb m = [2,3]. 

This equation is relevant because an order parameter for DCSB is the chiral- 
limit, vacuum quark condensate [2]: 



(qq)l := Z 4 (fi 2 ,A 2 ) N c J tr D [S Aocmbm=0 (g)] , (3) 



where Z 4 (/i 2 ,A 2 ) = [a(A 2 )/a:(/i 2 )] 7m(1+?/3) at one-loop order, with /i 2 the 
renormalisation point, £ the covariant-gauge fixing parameter [£ = spec- 
ifies Landau gauge] and 7 m = 12/(33 — 2Nf) the gauge-independent mass 
anomalous dimension. The ^-dependence of Z^fi 2 , M 2 ) is just that required 
to ensure that (qq)^ is gauge independent. It follows from (1) that an equiv- 
alent order parameter is X := B(p 2 = 0) . Chiral symmetry is dynamically 
broken when, with m^ m = 0, X ^ 0. 

Equation (2) is a nonlinear integral equation and the properties of its solu- 
tion depend on the kernel, which is constructed from D^ u (k) and T^(q,p). As 
summarised in Ref. [4], the connected, dressed-gluon 2-point function, D^^k), 
satisfies an oft analysed DSE. The qualitative conclusion of these DSE studies 
is that if the ghost-loop in the gluon DSE is unimportant, which is tautological 
in the ghostless axial gauges, then relative to the free gauge boson propaga- 
tor the dressed-gluon propagator is significantly enhanced in the vicinity of 
k 2 = 0, where it is a regularisation of l/k A as a distribution [5] . That en- 
hancement persists to k 2 ~ 1-2 GeV 2 , where a perturbative analysis becomes 
quantitatively reliable. 

The other term in the kernel of (1) is T^(q,p), the connected, irreducible, 
renormalised dressed-quark-gluon vertex: p, q are the momentum labels of the 
quark and antiquark lines, and the total momentum P : = p — q. The analogue 
of this vertex in QED has been much studied and it is argued [6] that it should 
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not exhibit particle-like singularities at P 2 = 0.Q The reasoning is simple: 
such singularities do not arise at low order in perturbation theory and hence 
such a vertex contradicts perturbation theory in any domain on which a weak 
coupling expansion is valid. 

Another arguably stronger reason is that a singularity of this type signals the 
existence of a massless bound state that mixes with the gauge boson, and such 
states have not been observed. This feature can be elucidated by considering 
the colour-singlet, axial-vector vertex in QCD, which is the solution of 



tu 



J tu n 



(4) 



where the matrix T H specifies the flavour structure of the vertex, x^aWi P) : = 
S(q + )Tf (q; P)S(q-) with q + := q + r/pP, g_ := q — (1 — rjp)P and P the 
total momentum, and S(q) := di&g[S u (q), Sd(q), S s (q), . . .]. In (4), K is the 
fully-amputated, quark-antiquark scattering kernel: by definition it does not 
contain quark-antiquark to single gauge-boson annihilation diagrams, such 
as would describe the leptonic decay of the pion, nor diagrams that become 
disconnected by cutting one quark and one antiquark line. 



In the chiral limit the solution of this equation is [3] 



H 



rf#;P) = ^-75 



ltl F R {k- P) + 7 • kk u G R (k; P) - <v K H R {k- P) 



(5) 



+ f*(*;P) + f H ^T 5 H (k;P) 



where: F R , G R , H R and fg, are regular as P 2 -> 0; P a f ^(k; P) ~ 0(P 2 ); and 
T 5 H (k;P) is the Bethe-Salpeter amplitude for a massless pseudoscalar bound 
state; i.e., T 5 H (k; P) satisfies the associated homogeneous Bethe-Salpeter equa- 
tion. The vertex is gauge covariant: the pole-position and which is the lep- 
tonic decay constant, are gauge invariant[2] and the bound state amplitude 
responds in a well-defined manner to a gauge transformation. In 3-flavour, 
massless QCD the poles in the axial-vector vertices correspond to the octet of 
Goldstone bosons. There should be no such singularities in the colour-singlet 
vector vertex, and this is verified in model studies [7] . 



A particle-like singularity is one of the form (P 2 )~ a , a G (0, 1]. In this case one 
can write a spectral decomposition for the vertex in which the spectral densities 
are non-negative. This is impossible if a > 1. a = 1 is the ideal case of an isolated, 
(5-function singularity in the spectral densities and hence an isolated, free-particle 
pole, a G (0, 1) corresponds to an accumulation, at the particle pole, of branch 
points associated with multiparticle production. 
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Similar observations apply to the fully-amputated, dressed-quark-gluon ver- 
tex, Tfj,[q,p). It satisfies an integral equation like (4) with the complication 
that, in addition to the term involving K, there are 3 other terms involving 
the scattering kernels for: q-q to 2-gluon, K 2g ; q-q to ghost- ant ighost, K ghgh \ 
and q-q to 3-gluon, K 3g . Recall that, by definition, none of these kernels contain 
single-gluon intermediate states. Hence a massless, particle-like singularity in 
this vertex signals the presence of a colour-octet bound state in one of the scat- 
tering matrices: M := K/[l - (SS)K); M 2g := K 2g / [1 - (D D) K 2g ]; etc. As no 
such coloured bound states have been observed, one must reject calculations 
or Ansatze for any of the Schwinger functions that entail a particle-like singu- 
larity in this vertex. The same objection applies to particle-like singularities 
in the fully-amputated, dressed-3-gluon vertex, and all like n-point functions. 
This anticipates the result of an estimate [8] of the 3-gluon vertex via a nu- 
merical simulation of lattice-QCD, which shows no evidence for a singularity 
of any kind. 

From a phenomenological perspective, a combination of D^ y {k) enhanced as 
described and without particle-like singularities is an excellent result, 

since it is sufficient to yield DCSB and confinement^ without fine-tuning [1]. 
It can also provide for a quantitatively accurate description of a wide range of 
hadronic observables [1,3,9], although this depends more on the detailed form 
of D^ v (k) and T^q.p). 

Does the phenomenology of the strong interaction require that the gluon prop- 
agator be strongly enhanced relative to the free gauge-boson propagator? In 
Landau gauge 

DM = U» - A(fc 2 ) , A(k 2 ) : = 1 V(k 2 ) , (6) 



and the question can be posed as: Do the observable phenomena necessarily 
require (Aq^ ~ 220 MeV) 

P(fc 2 )>>lforO<fc 2 <10A 2 ? (7) 



We do not have an answer but we can explore alternatives. The antithesis of 



2 Herein confinement means that the dressed-quark and -gluon 2-point functions do 
not have a Lehmann representation; i.e., do not have a spectral representation with 
a non-negative spectral density. This is a sufficient but not necessary condition. [10] 
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(7) is the extreme^ possibility that 

V{k 2 = 0) =0, V{k 2 ) <lVk 2 , (8) 

which was canvassed in Ref. [11]. Therein the ghost-loop contribution to the 
gluon DSE is retained and the Ansatze for the 3-gluon and quark-gluon vertices 
exhibit ideal particle-like poles [a = 1]. Since these poles are an essential 
element of the solution procedure then, in the absence of a physically sensible 
interpretation or explanation of them, one could simply reject this result. 

Alternatively, one can suppose that (8) is more robust than the procedure 
employed to motivate it and explore the phenomenological consequences of 
the conjecture [11]: Vs{k 2 ) := fc 4 /(fc 4 + 6 4 ), where b is a dynamically generated 
mass scale. Following this approach it was found that if there are no particle- 
like singularities in the quark-gluon vertex, T At (g,p), then Vs(k 2 ) is unable 
to confine quarks [12,13] and b must be fine-tuned to very small values [b < 
b c ~ Aqcd] if DCSB is to occur [12-14]. It is therefore apparent that (8) is 
phenomenologically difficult to maintain.^ 

Nevertheless, the hypothesis has been explored in studies [15] of the dressed- 
gluon 2-point function using numerical simulations of lattice-QCD. V(k 2 = 0) 
is necessarily finite in simulations on a finite lattice because of the inherent 
infrared cutoff. Thus one can only truly determine V(k 2 ~ 0) by considering 
the behaviour of the numerical result in both the countable limit of infinitely 
many lattice sites and the continuum limit. The form Vs{k 2 ) does not provide 
as good a fit to the lattice data as an alternative form, which in the countable 
limit is 

V ^k 2 ) := M2 - / fc2 (fc2fl2) , , < k 2 < 0.6/ a 2 ~ 50 A^ CD , (9) 



where 1/a ~ 2.0 GeV is the inverse lattice spacing, Z « 0.1, r) ~ 0.53, and 
M w 0.16 GeV. This takes the maximum value V L {k 2 = 21A^ CD ) = 13.6 and 
corresponds to a less extreme alternative to (7), which we shall characterise 



3 "Extreme" because it corresponds to a screening of the fermion-fermion interac- 
tion, as familiar in an electro dynamical plasma, rather than the antiscreening often 
discussed in zero-temperature chromodynamics. 

4 DCSB requiring b ~ is indicative of the dynamical evasion of (8) since Vs{k 2 ) — > 
1 rapidly for small values of b. We do not consider the possibility that an irreducible 
vertex has a non-particle-like singularity; i.e., a singularity of the form (k 2 )~ a , a > 1, 
as there is no indication of such behaviour in any study to date. 
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V(k 2 = 0) = 0, max(v{k 2 )) <C(10). (10) 



The feature V(k 2 = 0) = is critically dependent on whether M is nonzero, 
or not. It appears to be nonzero in the countable limit but, as emphasised in 
Ref. [15], the behaviour of M (and if) in the continuum limit is unknown. 

The phenomenological implications of (9) can be explored using the methods 
of Ref. [12]. A preliminary estimate follows by observing that "Pl(A; 2 ) is ap- 
proximately equivalent to Vs{k 2 ) if one identifies b^ ~ ^ M/a = 0.57 GeV. 
Hence one expects that (9) does not generate DCSB nor confine quarks.0 In 
order to quantitatively verify this conclusion we note that: it is the combi- 
nation g 2 V(k 2 )/k 2 that appears in (1) and g 2 is not determined in Ref. [15]; 
and we must extrapolate "Pl(/c 2 ) outside the fitted domain. Both of these re- 
quirements are fulfilled if we assume that a one-loop perturbative analysis is 
reliable for k 2 > 25 Aq CD , and employQ 

9 ^)J^ e) ' e ^ (u) 

requiring that Ai(k 2 ) := Vi(k 2 )/k 2 and its first derivative be continuous at 
k^. This procedure yields Ai(k 2 ) in Fig. 1 with 

g m = 0.Q5, k 2 m = 30 Aq CD . (12) 



It is now straightforward to solve (1) with a variety of Ansatze for the quark- 
gluon vertex that do not exhibit particle-like singularities.^ We employed the 
bare vertex T^p, q) := 7 M ; the Ansatz [6]: 

5 A dressed-gluon propagator satisfying (8) automatically satisfies (10). The model 
of Ref. [16] is in the class specified by (10), as are the fitted forms obtained in all 
existing lattice-QCD simulations. 

6 A value of b ~ 0.4GeV> b c in Vs(k 2 ) provides the best fit to the lattice data 
and this supports the same conclusion. 

7 In representing V(k 2 ) in the ultraviolet via the running coupling, g 2 (k 2 ), we are 
effectively enforcing the identity: Z\ = Z3 between the ghost-gluon-vertex and ghost 
wave function renormalisation constants. This "Abelian Approximation" is a phe- 
nomenologically well justified Ansatz that ensures [3] the correct one- loop anomalous 
dimension for both: the dressed-quark propagator obtained as a solution of (1); and 
the pseudoscalar meson bound state amplitudes. At most it introduces an error at 
large-A; 2 in the power of the ln-dependence of V(k 2 ). 

8 Equation (11) defines a renormalisable model quark DSE, which we solved in 
the manner described in Ref. [3]. For simplicity, we renormalised at the momentum 
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Fig. 1. Ai(k 2 ) := Vi{k 2 )/k 2 from (11). V\{k 2 ) is (9) in the infrared and extrapolates 
this lattice model outside the domain accessible in the simulation [15]. 



iff (p,g):=*£ A (p 2 ,g 2 )7 M 



(13) 



Hp + q)n \¥i -(p + q) & a (p 2 , g 2 ) + A B (p 2 , q 2 ) 



where E A (p 2 ,g 2 ) := [A{p 2 ) + A(q 2 )]/2, A A (p 2 ,q 2 ) ■= [A(p 2 



A(q 2 )]/[p 2 -q 2 ] 



and A B (p 2 , q 2 ) :— [B(p 2 ) — B(q 2 )]/[p 2 — q 2 ]; and an augmented form [17] 



rf(p;<f):=rf(p, 9 ) + rj(p, 9 ), (14) 
r»( P , ,) := l^fl^LtMLI^j} [m _ m] , (15) 

with d{p, q) := ([p 2 -g 2 ] 2 + [M(p 2 ) 2 + M(g 2 ) 2 ] 2 )/(p 2 + g 2 ), each of which allows 
the quark DSE to be solved in isolation. In all cases we found X = in the 
chiral limit; i.e., no DCSB with B(p 2 ) = 0. 

The absence of DCSB means it is straightforward to decide whether (11) gener- 
ates confinement. In this case quark confinement is manifest if Z(p 2 ) is smooth 
and vanishes at p 2 = 0; while the existence of a Lehmann representation and 
the concomitant lack of confinement is clear if Z(p 2 ) does not vanish at p 2 = 0. 
In Fig. 2 we plot the solution Z(p 2 ) obtained from (1) with the vertex Ansatze 
introduced above. The behaviour of the solution is qualitatively equivalent in 

cutoff, Auv ~ 10 4 Aqcd> since the p 2 -evolution of A(p 2 ) and B(p 2 ) beyond that 
point is completely determined by g 2 (k 2 ). 
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Fig. 2. Z(p 2 ) obtained as the solution of (1) using (11) with: (13) - solid line; (14) 
- dashed line; and r M (p;g) = 7^ - dotted line. That (9) does not confine quarks is 
manifest in the result: Z(p 2 = 0) ^ 0, which is independent of the vertex Ansatz. 

each case. This demonstrates explicitly that (9) neither generates DCSB nor 
confinement .0 thereby confirming our preliminary hypothesis based on the 
correspondence with Vs via an effective value of b. 

We have also solved the equation for Z{p 2 ) using 



Vi(k 2 



1 +^e 



-fc 2 /A, 



QCD 



V^k 2 



(16) 



where q is a variable "strength" parameter. Increasing q moves the peak in 
Vi(k 2 ) toward k 2 = and increases its height, thereby making it increasingly 
like the model of Ref. [3]. The form of Z(p 2 ) is qualitatively unchanged and 
hence there is no signal for the onset of confinement until q > 300. At s — 300 
the maximum value is 



Viik 



1 -0.98A| CD ) 



210 



(17) 



9 The result of a recent numerical simulation [18] is pointwise smaller in magnitude 
(< 1/3) than (9) on the entire fitted domain and hence the same conclusions apply 
in that case. The discrepancy is neither identified nor explained in Ref. [18]. Both 
simulations use (5 = 6.0. Ref. [15] has 500 configurations on a 24 3 x 48 lattice with 
(5 M A A1 (x))Lat. < 10~ 6 , while Ref. [18] has 75 configurations on a 32 3 x 64 lattice with 

(VW) Lat . < 10_12 - 
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and Vi(0.98A 2 QCD )/Vi(30A 2 QCD ) = 16, cf. Vi(21 A^ CD )/^(30 A| CD ) = 1.0. In 
the model of Ref. [3] the peak is at k 2 = 3.7 Aq CD and the value of this ratio is 
44, neglecting only for the purpose of this comparison the purely long-range, 
£ 4 (fc)-part of that interaction. A comparison of g 2 n T '«(Aq CD ) ~ 89 with the 
critical coupling of g 2 ~ 11 in Refs. [19] shows that such large values of <, 
ensure DC SB. 

The hypothesis (8) has also re-emerged recently in a DSE study [20] that is 
qualitatively akin to Ref. [11]: in its result that V(k 2 ) = Ak 4 for k 2 ~ 0; in 
postulating a significant role for the ghost-loop in the gluon propagator; and 
in employing a ghost-gluon vertex that is free of particle-like singularities!" 111 ! 
Following the above analysis, the results of Refs. [12,13] can be applied directly 
in this case. The gluon propagator is smooth and can be characterised by a 
value of 6 4 ~ A 4 A4, where A ~ 1 and A is a mass-scale that is left undetermined 
in Ref. [20]. Choosing any reasonable value of A; e.g., A > Aq CD , this gluon 
propagator, with a quark-gluon vertex that is free of particle-like singularities, 
neither yields DCSB nor confines quarks. 

It was recognised in Ref. [15] that the hypotheses (8) and (10) are problem- 
atic, and this is emphasised by our results: in the absence of particle-like 
singularities in the dressed-quark-gluon vertex, the possibility of a dressed- 
gluon propagator satisfying (10) is excluded by the existence of DCSB. In the 
lattice simulations the infrared scale, M, that entails (10) may vanish in the 
continuum limit: DCSB and confinement appear to require that. However, this 
limit is presently unexplored. In the DSE study the treatment of the ghost- 
gluon vertex and ghost-gluon scattering kernel is rudimentary and, with the 
information currently available to us, it is difficult to estimate the sensitivity 
of the results to the truncations. The incompatibility we identify between (8) 
and the phenomena of the strong interaction suggest that an examination of 
the effects of these truncations is necessary. 

Our primary predicate is that the dressed-quark-gluon vertex should not ex- 
hibit a particle-like singularity. If this is false and the vertex exhibits a singular- 
ity at P 2 = then (10), or even (8), may be reconcilable with the phenomena 
of the strong interaction. However, whether that is truly the case will likely 
depend on details and therefore require fine-tuning in the theory. 
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haviour is observed in lattice-QCD estimates of the 3-gluon vertex [8]. 
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